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II metodo comprende le fasi di: inviare in un volume sotto indagine una serie di 
segnali ultrasonici di eccitazione, tramite una cortina di trasduttori allineati lungo una 
direzione trasversale (jc), detti segnali ultrasonici propagandosi in profondita in detto 
volume secondo una direzione di propagazione (y); raccogliere, tramite detti trasduttori, 
segnali riflessi da riflettori disposti nel volume sotto indagine; eseguire su detti segnali 
riflessi una trasformata in direzione trasversale da un dominio spaziale (jc,y), definito da 
delta direzione trasversale (jc) e da delta direzione di propagazione (y), ad un primo 
dominio trasformato;. applicare, nel dominio trasformato, una trasformazione 
bidimensionale, per raddrizzare ogni immagine curva (Ipl; Ip2, Ip3) di un riflettore in detlo 
volume sotto indagine e renderla sostanzialmente ortogonale alia direzione di propagazione 
(y); comprimere, nella direzione trasversale (x) ciascuna curva raddrizzala (Iml, Im2, Im3) 
per concentrare detta immagine raddrizzata in una zona centrata in corrispondenza della 
posizione di detto riflettore lungo detta direzione trasversale (x). 
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aFirenze 

"Metodo di focalizzazione di una immagine ecografica e relativo sistema ecografico 

Descrizione 

Campo Tecnico 

La presente invenzione riguarda un metodo per la focalizzazione di immagini 
ecografiche, ed in particolare per eseguire la focalizzazione in senso traverso, cioe nella 
direzione parallela alio sviluppo della cortina di trasduttori formanti la sonda ecografica. 

L'invenzione si riferisce anche ad un sistema ecografico, cioe ad una 
apparecchiatura ecografica che attua detto metodo. 
Stato della Tecnica 

Nelle indagini ecografiche. si fa uso di sonde generalmente costituite da una serie di 
trasduttori od elementi piezoelettrici disposti su una cortina monodimensionale, oppure da 
un singolo trasduttore a scansione meccanica. Le ecografie sono prodotte affiancando le 
tracce provenienti da trasduttori diversi, o da posizioni diverse del singolo trasduttore nel 
caso di sonda monoelemento. I trasduttori generano segnali ultrasonici che si propagano nel 
corpo da sottoporre ad indagine e ricevbno il segnale ultrasonico riflesso dai riflettori che si 
trovano nel corpo stesso. Dal segnale riflesso o retropropagato si ricostruisce Pimmagine 
delFinterno del corpo esaminato. 

Per una ricostruzione perfetta della sezione indagata nell'ecografia, i trasduttori 
dovrebbero avere una sensibilita apprezzabile solo in un angolo solido infinitesimo, in 
modo da ricevere gli ultrasuoni provenienti da una sola direzione. In realta un trasduttore 
presenta un diagramma .di radiazione con un lobo di dimensioni non nulle attorno alia 
direzione principale ed eventualmente anche da lobi secondari. La Fig. 1 mostra un 
diagramma di radiazione di un trasduttore ultrasonico, con un lobo princ^ serie 

# 

Pag. 2 di 36 

C:\documenli_michele\descrizioni\acll91102.doc 




Uffici Tronic Ing. A. Mannucci srt 

di lobi secondari, rappresentati nello spazio tridimensionale x, y, z. 

Per questo motivo ogni trasduttore invia e riceve gli ultrasuoni non solo nella direzione 
principale (la direzione dell'asse z nella Fig.l) ma anche lateralmente. Nel piano di 
indagine della cortina, gli ultrasuoni retro-propagati da un riflettore presente nel tessuto 
sotto indagine sono awertiti non solo dal trasduttore piu vicino (cioe da quello che si trova 
allineato sull'asse z con il riflettore) ma anche da quelli circostanti, tanto piu numerosi 
quanto piu e grande il lobo di radiazione dei trasduttori. Nel segnale acquisito il riflettore 
compare pertanto in molte tracce e con ritardi diversi, proporzionali alia distanza tra il 
riflettore ed i vari trasduttori della cortina. In queste condizioni, neH'immagine costruita 
dall'ecografo, il riflettore non appare come un punto ma come una curva. 

Questo fenomeno e mostrato schematicamente in Figg.2 e 3. In Fig.2(a) si osserva 
come gli ultrasuoni emessi da un generico trasduttore Tr vengono propagati lungo direzioni 
diverse, mentre in Fig.2(b) si osserva come un generico trasduttore Tr riceve gli ultrasuoni 
riflessi da riflettori Rl, R2 disposti in posizioni diverse. In Fig.2(c), infine, e mostrato come 
il segnale eco prodotto da un singolo riflettore viene ricevuto non solo dal trasduttore Tr 
piu vicino ad esso, ma anche dai trasduttori adiacentL Quando un trasduttore riceve un 
segnale eco da un riflettore, tale segnale viene visto come se provenisse sempre da un 
riflettore disposto direttamente al di sotto del trasduttore stesso, cioe lungo l'asse di 
propagazione del lobo principale del trasduttore, anche quando si tratta in realta di un 
segnale propagato da un riflettore lateralmente spostato rispetto al trasduttore in questione. 
Inoltre, il trasduttore "vede" tale riflettore ad una profondita maggiore di quella reale. 

Poiche l'andamento delle distanze di un singolo riflettore rispetto ai vari trasduttori 
della cortina e iperbolico, Timmagine di tale riflettore che la cortina di trasduttori riceve e 
una curva costituita da un ramo di iperbole, come mostrato esemplificativamente nella 
Fig.3. Qui con S e indicata genericamente la sonda costituita da una cortina di trasduttori Tr 
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allineati lungo la direzione x, mentre con y e indicata la direzione di propagazione del 
fronte d'onda. II riflettore e localizzato nel punto di coordinate x 0 , yo, e viene "visto" dalla 
cortina di trasduttori come un ramo di iperbole C, con il vertice posizionato in 
corrispondenza del riflettore. 

E' evidente che le imrnagini che si ottengono in questo modo sono pesantemente 
alterate, non sono utilizzabili direttamente e devono essere elaborate. La focalizzazione in 
ecografia consiste esattamente in questo: riportare le imrnagini il piu possibile vicino a 
quelle che si sarebbero ottenute con trasduttori aventi lobo di radiazione di dimensioni 
infinitesime. 

Nel caso di sonde a cortina lineare, costituite da una cortina di trasduttori Tr disposti 
lungo una linea retta sono necessari due tipi di focalizzazione: in senso trasverso, cioe nella 
direzione (x) in cui si estende la cortina di trasduttori, e nella direzione perpendicolare ad 
essa ed alia direzione (y) di propagazione degli ultrasuoni, cioe nella direzione ortogonale 
al piano della Fig.3. La secohda focalizzazione si ottiene inserendo davanti alia cortina di 
trasduttori una lente acustica focalizzata airincirca a meta della profondita massima per cui 
si pensa di usare la sonda. 

Per quanto riguarda, viceversa, la focalizzazione in senso traverso, generalmente gli 
ecografi utilizzano sia una focalizzazione in trasmissione, come descritto di seguito, ma 
anche di una focalizzazione in ricezione, che e la focalizzazione vera e propria 
dePimmagine e che produce ogni traccia deH'immagine focalizzata sommando ad essa, 
dopo averle ritardate opportunamente, un certo numero di tracce ricevute da elementi 
adiacenti della cortina, come descritto piu avanti. 

La focalizzazione in trasmissione non e una focalizzazione deirimmagine, bensi si tratta 
di una tecnica per concentrare ad una certa profondita Fenergia trasmessa. Non viene usata 
una lente acustica ma ne viene simulato Teffetto elettronicamente. Ogni impulso 
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ultrasonico viene inviato quasi contemporaneamente da un gruppo di trasduttori adiacenti, 
rispetto a quello centrale. Ogni trasduttore del gruppo emette un impulso con un anticipo 
definito da una legge teoricamente iperbolica. Nella pratica, peraltro, vari ecografi 
utilizzano un'approssimazione parabolica 
5 In questo modo gli impulsi trasmessi si ritrovano in fase tra di loro ad una certa distanza 

dalla sonda; in questa posizione si dice che e stato posto il fuoco in trasmissione. 

Per concentrare l'energia trasmessa a varie distanze dalla cortina, si deve ripetere 
l'operazione descritta variando di volta in volta la posizione del fuoco. Questo 
procedimento pero richiede piii di un impulso per ottenere ogni singola traccia, quindi 
10 rallenta drasticamente la generazione delPimmagine completa. Per ottenere un'immagine 
ne vanno in effetti acquisite tante quanti sono i fuochi usati in trasmissione. Negli 
apparecchi commerciali generalmente e possibile posizionare da un minimo di un fuoco in 
trasmissione ad un massimo di pochi fuochi, per mantenere accettabile il numero di 
immagini al secondo presentate all'utente. 
15 La focalizzazione delle immagini ecografiche viene eseguita principalmente durante 

la ricezione degli ultrasuoni. Anche in questa fase non sono utilizzate lenti acustiche ma 
accorgimenti elettronici che ne simulano 1'effetto. Dopo la trasmissione di un impulso, gli 
ultrasuoni si propagano nel materiale sotto esame, vengono riflessi e tornano indietrb verso 
la sonda. Quelli che provengono da un riflettore posto ad una certa profondita raggiungono 
20 i trasduttori della cortina con tempi che seguono un andamento iperbolico. Anche gli 
impulsi elettrici generati da questi trasduttori sono ritardati tra loro con lo stesso 
andamento. La focalizzazione in ricezione consiste nel sommare al segnale elettrico 
prodotto da un trasduttore della cortina quelli prodotti da un gruppo di trasduttori adiacenti, 



dopo aver compensato i ritardi da cui sono affetti. Questo proc^to^to-viei^mpetuto per 
25 ogni traccia, spostando di un trasduttore il gruppo di trasduttor|gs^to^p^ ^g-J 
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Poiche la legge che regola i ritardi con cui gli ultrasuoni raggiungono la cortina 
varia con la profondita, anche quella usata per compensare tali ritardi deve variare a 
seconda della profondita che si vuol mettere a fuoco. Ma questo vuol dire che la curva dei 
ritardi usata nella focalizzazione in ricezione varia col passare del tempo, poiche nella fase 
di ricezione, via via che passa il tempo, arrivano alia cortina ultrasuoni provenienti da 
profondita sempre piu grandi. Per questo si parla di focalizzazione dinamica in ricezione. 
Tutto funziona come se davanti ad ogni trasduttore ci fosse una lente con focale variabile 
durante il tempo di ricezione. 

La corretta implementazione della focalizzazione dinamica in ricezione pone alcuni 
problemi pratici. Non e banale modificare con continuity la curva di ritardi usata per 
combinare le tracce, si pud pensare di discretizzarne le variazioni, suddividere la profondita 
esaminata in intervalli ed usare una sola curva in ciascun intervallo. I risultati ottenuti 
saranno tanto peggiori quanto tninore e il numero di curve diverse usate. In alcuni 
strumenti inoltre la curva dei ritardi non e un arco di iperbole ma un ramo di parabola che 
approssima la curva esatta: e questa un'ulteriore approssimazione del caso ideale. che 
peggiora i risultati. 

Ammettendo di riuscire ad implementarla in modo corretto, se Tunica 
focalizzazione utilizzata fosse quella in ricezione, le immagini ecografiche potrebbero 
essere elaborate ottenendo buoni risultati. Tuttavia, generalmente gli ecografi affiancano 
anche la focalizzazione in trasmissione, per far fronte alTattenuazione che le onde 
acustiche subiscono propagandosi nei tessuti organici. In queste condizioni, gli ultrasuoni 
arrivano ai trasduttori con la morfologia prevista dalla focalizzazione in ricezione solo se 
prdvengono dalle profondita corrispondenti ai fuochi in trasmissione. Gli ultrasuoni 
provenienti da altre profondita sono alterati dalla focalizzazione in trasmissione e non 
possono essere correttamente elaborati. In questo modo si ottiene una risoluzione non 
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uniforme con la profondita, buona solo in corrispondenza dei fuochi usati in trasmissione. 

Una focalizzazione non accurata delle ecografie, oltre ad una limitata risoluzione, 
puo generare anche artefatti nelle immagini, owero dettagli che in realta non sono presenti 
nella sezione indagata. 

L'effetto piu visibile e comunque Tallargamento di ogni punto dell'immagine in 
una figura che ha forma edimensioni dipendenti dalla profondita. 
Scopi e sommario deirinvenzione 

Scopo della presente invenzione e la realizzazione di un metodo per la 
focalizzazione in senso traverso, cioe parallelamente airallineamento liingo cui sono 
disposti i trasduttori formanti la cortina di trasduttori della sonda, di una immagine 
ecografiea che non presenti gli inconvenienti delle tecniche tradizionali. 

Piu in particolare, scopo della presente invenzione e la realizzazione di un metodo 
di focalizzazione e di un sistema ecograf ico che utilizza tale metodo, che raggiunga risultati 
qualitativamente migliori rispetto ai metodi conosciuti. * 

Sostanzialmente, il metodo secondo l'invenzione prevede le seguenti fasi: 

• inviare in un volume sotto indagine una serie di segnali ultrasonici, tramite una 
cortina di trasduttori allineati lungo una direzione trasversale, i quali segnali 
ultrasonici si propagano in profondita nel volume secondo una direzione di 
propagazione; 

• raccogliere, tramite detti trasduttori, i segnali riflessi da riflettori presenti nel 
volume sotto indagine; 

• previa demodulazione coerente del segnale per riportarlo in banda base, owero per 
traslarne lo spettro a cavallo della frequenza zero, eseguire sui segnali riflessi una 
trasformata in direzione trasversale da un dominio spaziale, definitp dalla direzione 
trasversale e dalla direzione di propagazione, ad un primo dominio trasformato; 
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• applicare, nel dominio trasformato, una trasformazione bidimensionale, per 
raddrizzare ogni immagine curva di un riflettore e renderla sostanzialmente 
ortogonale alia direzione di propagazione; 

• comprimere, nella direzione trasversale ciascuna curva raddrizzata per concentrare 
P immagine raddrizzata in una zona centrata in corrispondenza della posizione di 
detto riflettore lungo la direzione trasversale. 

Questa ultima fase comprende, in pratica, anche una antitrasformata di Fourier che 
riporta 1' immagine nel dominio spaziale 

In pratica, il metodo secondo l'invenzione si basa sulla compressione d'impulso in 
direzione laterale del segnale acquisito con una sbrida che non utilizza la focalizzazione 
elettronica, operazione preceduta da una demodulazione coerente dei segnali acquisiti. 

Per poter operare la compressione d'impulso in direzione laterale, essendo le 
immagini eeografiche affette da una pesante curvatura come esposto in cio che precede, 
F immagine e stata portata nel dominio della frequenza tramite una trasformata nella 
direzione ortogonale a quella di propagazione deirimpulso. Nel dominio trasformato si 
esegue una rimappatura che "raddrizza" V immagine spostando i suoi punti lungo la 
direzione parallela a quella di propagazione deirimpulso nel volume sotto esame. Tramite 
la rimappatura viene completamente eliminata la curvatura dell' immagine e diviene quindi 
possibile comprimere Timpulso in senso laterale, ottenendo cosi un'immagine focalizzata. 

In una pratica forma di realizzazione, la trasformata nel senso trasversale pud essere una 
trasformata di Fourier. In considerazione della complessita della funzione che esprime il 
segnale delPimmagine, tale trasformata puo essere in realta approssimata con il principio 
delle fasi stazionarie. La compressione del segnale in direzione trasversale consiste nei tre 
passaggi: trasformata, compensazione della fase e del gudagno, antitrasformata. La 
rimappatura consente di operare questa compressione deirimpulso tramite elaborazioni 
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eseguite su linee ortogonali alia direzione di propagazione degli ultrasuoni. 

II segnale ultrasonico che i trasduttori inviano, sequenzialmente, nel volume sotto 
indagine pud essere un segnale ecografico tradizionale. Alternativamente, e 
preferibilmente, peraltro, vengono usati segnali chirp, o comunque in generale segnali 
5 mpdulati in frequenza. Questo consente di raggiungere una maggiore penetrazione grazie 
alia maggiore energia del segnale ed una migliore definizione deU'immagine. 

Quando viene utilizzato un segnale modulato in frequenza, ad esempio con una 
modulazione lineare, occprre eseguire una compressione del segnale nel verso della 
propagazione, oltre che in direzione trasversale. Tale compressione viene 
10 vantaggiosamente - anche se non obbligatoriamente - eseguita con una trasformata tempo- 
frequenza. Anche in questo caso puo essere usata una trasformata di Fourier. 
Preliminarmente, peraltro, i segnali riflessi dai riflettori e captati dai trasduttori della sonda 
vengono riportati in banda base tramite una demodulazione coerente. 

La compressione del segnale nella direzione di propagazione, inoltre, comprende 
15 una operazione di compensazione della fase e del guadagno del segnale stesso. 

Come risultera chiaro dall'esposizione analitica del metodo oggetto della presente 
invenzione, Timmagine del volume sotto indagine viene costriiita inviando il segnale 
ultrasonico dai vari trasduttori che vengono eccitati in successione. In pratica, quindi, i vari 
trasduttori della cortina formante la sonda vegono eccitati l'uno dopo l'altro, ed 
20 analogamente accade in ricezione. Un trasduttore viene eccitato, trasmette un segnale 
ultrasonico e si mette in condizioni di ricezione del segnale ecografico. Dopo un tempo 
fissato, viene eccitato il trasduttore adiacente lungo la cortina e cosi via. Cio differenzia il 
metodo descritto rispetto alia focalizzazione in ricezione degli ecografi tradizionali, che 
utilizza molti trasduttori per ricevere il segnale riflesso e somma isegnali elettrici ottenuti 
25 per costruire una traccia deU'immagine. 
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SullMnsieme dei segnali cosi raccolti, che rappreseritano l'immagine iniziale, 
vengono eseguite le operazioni di focalizzazione. 

L'invenzione riguarda, anche un sistema ecografico comprendente almeno una 
sonda con una cortina di trasduttori allineati secondo una direzione trasversale di 
5 allineamento, mezzi per eccitare sequenzialmente detti trasduttori, mezzi di elaborazione 
per ricevere ed elaborare segnali riflessi da riflettori contenuti in un volume sottq indagine 
in cui vengono propagati segnali ultrasonici generati da detti trasduttori, detti mezzi di 
elaborazione eseguendo la focalizzazione in direzione trasversale del segnale ecografico 
ricevuto da detti trasduttori con un metodo secondo una o piu delle rivendicazioni 
10 precedenti. 

Breve descrizione dei disegni 

II trovato verra meglio compreso seguendo la descrizione e Tunitb disegno, il quale 
mostra una pratica forma di attuazione non limitativa dell'invenzione. Piu in particolare, 
nel disegno mostrano: la . 
15 Fig.l un diagramma di radiazione di un trasduttore ultrasonico; la 

Fig.2 una rappresentazione schematica dell 5 effetto della sensibilita laterale di un 
trasduttore di una cortina lineare di trasduttori di una sonda ultrasonica; la 

Fig.3 uno schema che mostra una sonda ultrasonica e rimmagine di un riflettore 
puntiforme da essa rilevato, in assenza di focalizzazione; la 
20 Fig.4 un diagramma mostrante Tandamento indicativo della sensibilita laterale di un 

singolo trasduttore di una cortina di trasduttori; la 

Fig.5 l'andamento indicativo dell'inviluppo di un segnale modulato linearmente in 
frequenza, emesso da uh trasduttore; la 

Fig.6 una rappresentazione a blocchi di un demodulatore coerente; la 
25 Fig.7 una rappresentazione schematica di un segnale generato da un riflettore 
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singolo, prima della focalizzazione; la 

Fig.8 una rappresentazione schematica dei segnali generati da tre riflettori nello 
spazio; la 

Fig. 9 una rappresentazione dei segnali generati dai tre riflettori di Fig.8 nel dominio 
5 trasformato; la 

Fig.10 il segnale generato dai riflettori di Fig.8, nel dominio trasformato, dopo la 
rimappatura; la 

Fig. 11 il segnale generato dai tre riflettori di Fig.8 dopo Tantitrasformata in senso 
traverso; la 

10 Fig.12 uno schema a blocchi riassuntivo del metodo secondo l'invenzione in una 

forma di realizzazione che fa uso di un segnale modulato in frequenza; la 

Fig.13 lo stesso schema nel caso di segnali non modulati in frequenza; la 
Fig. 14 i segnali chirp usati per le prove sperimentali; la 
Fig.15 i segnali non modulati usati per le prove sperimentali; la 

15 Fig.l6(a) una immagine reale acquisita dall'ecografo di un filo teso in acqua; la 

. Fig:16(b) la traccia centrale deirimmagine di Fig. 16(a); la 

Fig.l7(a) la parte reale deirimmagine di Fig.l6(a) demodulata coerentemente; la 
Fig.l7(b) ringrandimento della traccia centrale deirimmagine di Fig.l7(a); la 
Fig. 18(a) il modulo deirimmagine di Fig. 17(a) focalizzata in profondita; la 

20 Fig.l8(b) un ingrandimento della zona centrale dell'immagine di Fig.l8(a); la 

Fig.l9(a) Tandamento del modulo della traccia evidenziata hell' immagine di 
Fig.l8(a), normalizzato ad 1; la 

Fig.l9(b) la stessa traccia di Fig.l9(a) rappresentata in scala logaritmica; la 
Fig.20(a) il modulo della trasformata di Fourier in senso traverso deirimmagine 

25 focalizzata in profondita di Fig. 18(a); la 
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Fig.20(b) un ingrandimento di una porzione della Fig.20(a); la 
Fig.21(a) Timmagine di Fig.20 dopo la rimappatura; la 
. Fig.21(b) un ingrandimento di un dettaglio centrale di Fig.21(a); la 
Fig.22(a) il modulo deirimmagine focalizzata; la 

Fig.22(b) la stessa immagine di Fig.22(a) rappresentata con una scala di grigi 
logaritmica; la 

Fig.23(a) il modulo della traccia orizzontale, centrata sul filo, deirimmagine 
focalizzata, normalizzato ad 1; la 

Fig.23(b) lo stesso modulo di Fig.23(a) in scala logaritmica; la 

Fig.24 il segnale reale acquisito dalPecografo, generato da due fili tesi in acqua, uno 
in metallo e 1'altro in nylon; la 

Fig.25 il modulo deirimmagine di Fig.24 dopo la focalizzazione, con un 
ingrandimento di dettaglio; la 

Fig.26 il segnale reale acquisito dall'ecografo generato da un filo teso in acqua, 
ottenuto con un segnale non modulato in frequenza; la 

Fig.27(a) il modulo deirimmagine di Fig.26 dopo la focalizzazione, in scala liheare 
di grigi; la 

Fig.27(b) lo stesso modulo deirimmagine di Fig.26 dopo la focalizzazione, in scala 
logaritmica dei grigi; la 

Fig.28(a) il modulo della traccia orizzontale, centrata sul filo, deirimmagine 
focalizzata di Fig.27, normalizzato a 1; la 

Fig.28(b) lo stesso modulo di Fig.28(a) in scala logaritmica; la 

Fig.29(a) il modulo della traccia verticale centrata sul filo deirimmagine focalizzata 
di Fig.27, normalizzato a 1; la 

Fig.29(b) lo stesso modulo di Fig.29(a) in scala logaritmica; le 
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Figg.30(a) e (b) una immagine focalizzata tradizionalmente di un filo teso in acqua, 
in una scala lineare di grigi ed in scala logaritmica, rispettivamente; le 

Figg.31(a), (b) e (c) il confronto del modulo delle tracce orizzontali centrate sul filo 
ottenuto con una focalizzazione tradizionale e con la focalizzazione del metodo della 
5 presente invenzione, senza e con modulazione del segnale, rispettivamente; e le 

Figg.32(a), (b) e (c) lo stesso confronto di Figg. 31 sulla traccia verticale centrata 
sul filo. 

Descrizione dettafiliata delle forme di attuazione preferite deir invenzione 
1. Descrizione del metodo di elaborazione 
10 In Fig.3 e illustrata schematicamente una sonda ecografica S con una cortina di 

trasduttori Tr allineati lungo una direzione x, detta nel seguito anche direzione principale 
della cortina di trasduttori. Nel metodo della presente invenzione i trasduttori Tr vengono 
eccitati in successione per emettere un segnale ecografico. Nell'esempio qui illustrate, il 
segnale emesso e modulato in frequenza, cioe in pratica e un segnale cosiddetto chirp. II 
15 segnale riflesso dai riflettori contenuti nel tessuto - o comunque nel corpo sotto indagine - 
per ogni segnale ultrasonico emesso da un singolo trasduttore incide su molti trasduttori 
della cortina, che lo awertono con differenti ampiezze in base alPangolo di incidenza, ma 
per generare Timmagine viene registrato solo il segnale elettrico prodotto dal trasduttore 
che ha trasmesso Timpulso. 
20 La sensibilita laterale dei singoli trasduttori formanti la cortina ha un andamento a 

campana con un mjassimo in corrispondenza del proprio asse ed un andamento decrescente 
fino a zero lungo la direzione trasversale, come rappresentato in Fig A 

Un riflettore PI di dimensioni puntiformi, cioe infinitesime^chesi trova nel volume 
sotto indagine e che viene investito dal fronte d'onda generat^td^un sii%owtosduttore Tr 



25 posto in posizione x 0 genera una eco che viene ricevuta dai|singoli) trasduttof i-^ella cortina 
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lineare che vedono il riflettore PI entro un angolo O ma x. 

II segnale ricevuto da trasduttori che vedono il riflettore PI sotto un angolo 
maggiore di O max si suppone nullo, cioe la sensibilita si assume nulla al di fuori di tale 
angolo. Supponendo di lavorare con quantita continue, non ancora campionate, rimmagine 
che il riflettore PI genera sulla cortina di trasduttori Tr e analiticamente esprimibile come 

img(x,t ; x 0 , r) = a(x ; x Q , r )inv^t- ~i?(x;x 0 ,r)j- 

. C os|2^/ T +^^ (1) 

dove: 

x e la coordinata relativa alia direzione principale della cortina, cioe la direzione lungo cui 
si sviluppa la cortina di trasduttori; 
t e il tempo 

x 0 e / indiviuuaiio la posizionc del riflettore PI in senso traverso (lungo. x) ed in profondita 
(lungo la coordinata y in Fig.4) 

a(x,x 0 ,r) e la funzione di sensibilita laterale di ogni singolo trasduttore Tr della sonda. 
Questa sensibilita e rappresentata in Fig.4 dalla curva a campana rappresentante la 
sensibilita del trasduttore disposto in posizione x 0 . Essa e funzione di x 0 , r e x in quanto la 
sensibilita del trasduttore in posizione generica x al segnale riflesso dal riflettore PI 
dipende dalla posizione x 0 in cui si trova il trasduttore che ha emesso il segnale ultrasonico 
di eccitazione e dalla profondita a cui si trova il riflettore PI, profondita che e funzione 
diretta del tempo t intercorrente tra Tistante di emissione del segnale di eccitazione e 
Tistante di ricevimento del segnale riflesso da parte del generico trasduttore in posizione x; 
mv(t) e Tinviluppo deH'impulsd ultrasonico trasmesso. A seconda dei sistemi utilizzati 
potra" essere una finestra a campana o rettangolare, in ogni caso viene assunto nullo 
ovunque tranne che nell'intervallo [ -T/2, 772 ]. In Fig.5 e rappresentato un esempio di 
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forma d'onda del segnale di eccitazione, che pud essere utilizzato. In questo caso 
Tinviluppo ha forma rettangolare; 

c e la velocita dell'onda ultrasonica nel mezzo di propagazione, in specie nel tessuto sotto 
indagine; 

5 R e la distanza effettiva tra il generico trasduttore in posizione x ed il riflettore PI. Essa e 
data da 



R{x;x 0 r)=^(x-x 0 f + r 2 

(f>o e la fase iniziale del segnale, indipendente da x\ 

T e la durata deU'impulso ultrasonico emesso dal trasduttore Tr. Tipicamente tale durata e 
10 dell'ordine di 20 microsecondi. Questa durata, relatiyamente elevata, del segnale e dovuta 
al fatto che nell'esempio qui esposto (come sopra accennato) il segnale utilizzato e un 
segnale modulato in frequenza. Negli ecografi tradizionali, che usano (come verra 
accennato in seguito) segnali non modulati in frequenza, il segnale ha durata molto 
inferiore, dell'ordine del microsecondo. 
15 f T e la frequenza centrale del segnale di eccitazione, cioe del segnale ultrasonico emesso dal 
trasduttore in posizione x 0 ; 

A/e il range ehtro cui varia la frequenza / del segnale di eccitazione emesso dal trasduttore 
in posizione xq. Pertanto, neH'intervallo temporale T durante il quale viene emesso il 

A -f 

segnale ultrasonico, la frequenza del segnale stesso varia da una frequenza iniziale f T - — 

A -f 

20 ad una frequenza finale f T + — . 

- & 

Per maggior chiarezza Tespressione (1) del segnale rappresentante Timmagine 
,\ ; \ 1 [ ricevuto dalla cortina di trasduttori puo essere riscritto come: 

.' . J »" 

? iimg(x,t;x 0 ,r) = a(x,x 0 ,r)inv\t R(x;x 0 ,r)Y 

■ \ 0 ) 
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•cos 2;r/ T ^-^/j(x;Xo,^^ + p 0 



(2) 



A/" 

avendo posto b = n — . 

II segnale descritto dairespressione (2) e un'onda sferica, generata dal riflettore 
puntiforme PI ed incidente sulla cortina di trasduttori formanti la sonda, che si propaga a 
5 meta della velocita effettiva degli ultrasuoni nel mezzo sotto indagine: trattare questo 
segnale equivale a considerare la reale situazione del segnale che viene trasmesso dal 
trasduttore Tr* incide sul riflettore PI e ritorna indietro a velocita c. 

Per poter essere elaborato, il segnale ricevuto deve essere riportato in banda base 
tramite una demodulazione coerente. Cio yiene ottenuto con un filtro rappresentato in 
10 Fig.6. In sostanza con questa operazione di demodulazione si trasla la frequenza centrale 
sullo zero e mantenendo Tinformazione relativa alia fase. II segnale che si ottiene e 
espresso dai 

V C > ; . (3) 

NeU'espressione (3) vengono usati i termini gia definiti con riferimento alle 
15 precedent equazioni (1) e (2), ed inoltre e stata introdotta una costante di propagazione k, 

' 2n f 
p ar i a — owero pari a 2n— , essendo X la lunghezza d'onda del segnale ultrasomco, ed 

A, c 
essendo - 
1 



^ demod ^ ^ 



M 



NeU'espressione (3) si notano due contributi di fase rappresentati dai due termini 



20 ^ 'esponenziali: il primo e dovuto alia geometria del problema e dipende solo dalla coordinata 



x,il secondo alia modulazione di frequenza del segnale trasmesso e dipende da entrambe le 
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coordinate, xtt. 

Ora che il segnale e demodulato, conviene cambiare la coordinata temporale t con 
quella spaziale y che rappresenta la distanza nel verso di propagazione del fronte d'onda 

2 

ultrasonico, tenendo conto che f = — y> e quindi passare dalla funzione inv(t) a 

c 

5 p(y) « ' nv ^ t j e dalla costante b a b = -^-b . Con questo cambio di coordinata il segnale 
e definito dall'espressione: 

^,y;w) = ^«(*; W ' (4) 

II passaggio dalla coordinata temporale alia coordinata spaziale e potrebbe essere 

effettuato anche in una fase precedente deU'elaborazione. Peraltro, eseguirla dopo la 
10 demodulazione e vantaggioso ai fini della precisione del calcolo. Infatti, mentre. la 

coordinata temporale viene misurata dallo strumento, la coordinata spaziale e ottenuta da 

quella temporale tramite il legame spazio-temporale rappresentato dalla velocita di 

propagazione c degli ultrasuoni nel mezzo sotto indagine. Tale velocita non e conosciuta in 

modo esatto, bensi solo stimata e pud variare ad esempio in funzione del tipo di tessuto in 
15 cui l'onda ultrasonica si propaga. E' quindi opportuno eseguire la demodulazione, che 

utilizza le funzioni trigonometriche seno e coseno nella variabile temporale /, prima del 

cambio di coordinate. 

Poiche il segnale generato dal trasduttore Tr e un segnale modulato in frequenza di 

durata T relativamente elevata, delFordine di 10-20 microsecondi e quindi non ha la forma 
20 di un impulso, il segnale (4) che, riflesso dal riflettore PI raggiunge la cortina di trasduttori 

della sonda S ha una elevata durata e quindi, in assenza di adeguate elaborazioni, il 
, -<f»-triflettore PI verrebbe visto, su ciascuna traccia, come una striscia di lunghezza nel verso 
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corrispondente alia coordinata temporale t. La compressione dell'impulso comprende 
anche una operazione di compensazione della fase e del guadagno 

Benche sia possibile eseguire le operazioni suddette di compensazione della fase e 
del guadagno e di compressione in profondita nel segnale nel dominio del tempo, il calcolo 
risulta estremamente semplificato se eseguito nel dominio delle frequenze. A tale scopo il 
segnale demodulate, espresso dalla equazione (4), viene sottoposto ad una trasformata di 
Fourier in profondita, cioe nella coordinata y. L'operazione da eseguire e espressa da: 

pP(x 9 a> y ;x 09 r}*jl2 ... - ^ 

Data la complessita della funzione pp(x, y; xo , r) (definita dall'espressione (4)) 
1' integrate non e gestibile in forma esplicita. E' tuttavia possibile approssimarne il valore 
applicando il Principio della Fase Stazionaria o Metodo della Fase Stazionaria. Questo 
principio, conosciuto anche come Metodo della Fase Stazionaria, e una tecnica per 
approssimare integrali di funzioni complesse aventi fase che varia molto piu rapidamente 
del modulo. Supponendo di avere la funzione f(x), di cui vengono messi in evidenza il 
modulo g(x) e la fase e di volerne calcolare l'integrale, I: 

f(x)=g(x)e»V • 
/ = + Jg(x)e jHx) (Lc 

-cc • _ 

si pud dimostrare che tale integrale e approssimato daH'espressione: 



dove 

d0(x) 



"S/x tale che 

dx 
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1 B 2 *(x) 



*0. 



2 dx 2 

Per una trattazione analitica del principio si rimanda alia seguente letteratura: 

F. M.Henderson e A.J.Lewis, Principles and Applications of Imaging Radar , John Wiley & 
Sons, New York, 1998 

5 J.C.Curlander e R.N.McDonough, Synthetic Aperture Radar: Systems and Signal 
Processing, John Wiley & Sons, New York, 1991 

G. S.Kino, Acoustic Waves: Devices, imaging, and analog signal processing, Prentice-Hall, 
Englewood Cliffs, New Jersey, 1987 , 

M.Born e E.Wolf, Principles of Optics 6th ed., Pergamon Press, Oxford, 1980 (appendice 
10 III) 

E.ChassanderMottin e P.Flandrin, "On the stationary phase approximation of chirp spectra 
", in Proceedings of the IEEE-SP International Symposium on Time-Frequency and Time- 
Scale Analysis, Pittsburgh, PA, USA, pp. 117-126, 6-9 Ottobre 1998 

E-.Chassande-Mottin e P.Flandrin, "On the Time-Frequency Detection of Chirps", Applied 
15 and Computational Harmonic Analysis, vol. 6, n. 2, pp. 252-281, Marzo 1999 
A.Papoulis, Signal Analysis, McGraw-Hill, New York, 1977. 

Applicando dunque il principio della fase stazionaria al problema in 
discussione, il calcolo fornisce il seguente risultato: 



(6) 



I n jsegno(b)?- 

20 con C PFSy = J— e 4 



4 

"v, I tre termini esponenziali che compaiono nell'espressione (6) sono i termini di fase. 
II termine di fase lineare contiene l'informazione relativa alia posizione (nel dominio delle 

V\ . . . 
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frequenze) del riflettore, mentre il termine quadratico, che deriva dalla modulazione di 
frequenza del segnale trasmesso dal trasduttore Tr, deve essere eliminato dalla successiva 
operazione di compressione in profondita. In questo modo viene ottenuto lo spettro di un 
impulso molto stretto nella direzione y, posto alia profondita R(x; x 0 , r). Nella stesso 
passaggio e possibile compensare anche il fattore complesso introdotto dalla costante 
CpFsy e dalla prossima antitrasformata di Fourier. . 



e 4 J — 

Moltiplicando, quindi, il segnale espresso dalla (6) per la quantita , e 4b 

. 2y]7r\b\ 

si ottiene la nuova espressione: 



pP c (x,a> y ;x 03 r)=^C rfemo ^(x;x 0 ,r)/?f^-^ 



e - j2kR(x;x 0 ,r) ^- j a> y R(x,x a ,r) ^ 



Questa funzione rappresenta il segnale compresso in profondita, cioe nella direzione 
y. Infine, viene eseguita una operazione di antitrasformata di Fourier in profondita. Tale 
antitrasformata e facilmente calcolabile, ora che compare una dipendenza da co y 
semplicemente lineare sia per la fase che per il modulo, mentre e stato eliminato il termine 
quadratico. Eseguento 1' antitrasformata di (7) si ottiene: 

PPcfc^w)^^ (8) 
dove con p IFT (w) viene indicata la trasformata inversa di Fourier della funzione p(y). 
L' espressione (8) puo essere riscritta tramite la sostituzione p lFT (w)= p' lFT (2bw : 

PP c (x,y;x 0 ,r)=C de ^a(x;x 09 r)p^(y (9) 
Come p(y), anche p !FT (y) (antitrasformata di p(y)) e una finestra centrata sullo zero, ma 
di dimensione molto piu piccola. Poiche p /FT (y) si annulla per tutti i valori di y che non 
sona prossimi a zero, essa si awicina ad una delta di Dirac e l'espressione (9) .e 

sostknzialmente nulla ovunque tranne che suH'iperbole y = R(x;x Qy r) pesata dalla funzione 

• v 
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a(x; xo, r). Ne risulta un'immagine formata da un'unica curva, come si puo vedere in Fig.7. 

In sostanza, quindi, Poperazione di compressione in profondita, eseguita nel 
dominio delle frequenze previa trasformata di Fourier in y, comprime il segnale ricevuto da 
ciascun trasduttore della cortina formante la sonda S come segnale riflesso dal singolo 
5 riflettore PI, in un singolo punto dell'iperbole y - R(x;x 09 r). 

II segnale in questa forma non e ancora utilizzabile per ricostruire Pimmagine del 
volume sotto indagine. Infatti, in esso il singolo riflettore puntiforme PI e rappresentato da 
una iperbole con il vertice nella . posizione in cui si trova il riflettore stesso. Occorre 
eseguire una compressione in senso trasverso del segnale, owero nella direzione parallela 

10 alPasse x, cioe parallela alia cortina di trasduttori della sonda S. Tuttavia, poiche ciascun 
riflettore genera una immagine costituita da un ramo di iperbole caratterizzata da una 
marcata curvatura, se la compressione in senso trasverso avvenisse senza ulteriori 
elaborazioni preliminari, ne risulterebbe una perdita di definizione in profondita. Inoltre, 
nel dominio dello spazio i rami di iperbole costituiti dall' immagine generata da piu 

15 riflettori puntiformi a diverse profondita presentano curvature crescenti al diminuire della 
distanza dalla cortina di trasduttori, cioe curvature tanto piu marcate quanto minore e la 
coordinata r a cui si trova il riflettore. Questo significa che le varie iperboli generate da piu 
riflettori si intersecano con una cbnseguente impossibility di distinguere i contributi dei 
diversi riflettori. Questo fenomeno e rappresentato schematicamente in Fig. 8, dove tre 

20 riflettori puntiformi PI, P2, P3 sono disposti nel piano di giacitura degli assi di 
propagazione dei lobi di emissione principali dei trasduttori della cortina che si sviluppa, 
lungo la direzione x. Con n, x 0 r, r 2 , x 0 2 ed r 3 , x 03 sono indicate le coordinate dei tre 
riflettori. Con Ipl, Ip2 ed Ip3 sono indicati i tre rami di iperbole generati dai segnali riflessi 
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quanto si dispongono lungo Passe y secondo una disposizione con una curvatura crescente 
con y. La Fig.9 mostra il segnale generato dagli stessi tre riflettori nel dominio trasformato 
(&>x\ y). Con Ipl, Ip2, Ip3 sono ancora indicate le iperboli generate dai tre riflettori posti nei 
punti di coordinate v u x 0 i; r 2 , x 02 ed r 3 , x 0 3 nel dominio dello spazio (x, y). 

Pertanto, operazione preliminare alia compressione in senso trasverso del segnale e 
la trasformata di Fourier in senso trasverso del segnale compresso in profondita 
rappresentato dalPespressione (9): 

Come per la trasformata in profondita, anche per eseguire la trasformata di Fourier in 
senso trasverso si pud fare ricorso al Principio della Fase Stazionaria da cui si ottiene: 

Pp c («>, , y i *o ,r) - C PFS x C iemod a^ - ' ° > ' Rm ^' ,r) ;r y„ (y - R. {co x ; r)\ 



avendo definito la nuova funzione R^co^r) = . r (H) 



f- 



o 1 

X 

4k 2 



e la quantita C 



PFS x 



. (12) 



I CO \ 

2k\l- 



Ak 2 ) * 

Anche nelFespressione (10) della trasformata compare Pinviluppo p IFT (y) del 
segnale compresso in profondita, che e non nullo solo in prossimita di zero. Pertanto si puo 
con buona approssimazione affermare che nel dominio trasformato (o^; y) il segnale e 
apprezzabile solo sulla curva y = R„(co x \r). Tale curva e centrata sull'asse delle pulsazioni 
(co) \ dipende, oltre che da co x , unicamente dal valore di r. 

1 La curvatura di ciascun ramo di iperbole nel dominio trasformato (co*; y) significa 
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che il contributo di ciascun riflettore non si trova ad un'unica profondita lungo Passe y, 
bensi su una curva la cui esterisione in profondita non e trascurabile. Non e, quindi, 
possibile effettuare direttamente la compensazione di fase per la compressione in senso 
trasverso. Occorre prima una ulteriore elaborazione dell'immagine. 
5 Come gia anticipate, se sono presenti piii riflettori, a diverse posizioni in profondita 

(cioe con coordinate y diverse) e in senso trasverso (cioe con coordinate co* diverse), 
ognuno di questi da origine ad una curva che, a causa dell'espressione di R Q) {co x ;r) > ha 
curvatura diversa dalle altre ed e da esse distinta, senza sovrapposizioni, in quanto la 
curvatura delle iperboli cresce con y (Fig.9). In altri termini, per ogni punto del dominio . 
10 trasformato {cc^\ y) passa una sola curva qualunque sia il numero di riflettori puntuali che le 
hanno generate. 

E' allora possibile applicare una trasformazione airimmagine definita 
dall'espressione (11) per riportare ogni curva, caratterizzata da un certo valore di r, ad una 
linea retta ady costante (cioe parallela all'asse co* delle pulsazioni) e uguale al valore di r. 
15 Tramite l'equazione y - 'R a (a) x ;r) e possibile risalire a r conoscendo co x e y: 



(13) 

Si defnisce la trasformazione bidimensionale 

5) x =<o x 

(14) 



20 



che applicata all'immagine (10) produce: 



y-r 



L'espressione (15) e semplificabile, senza commettere errori significativi grazie 
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all'elevato valore di k, considerando le seguenti approssimazioni: 



Da cui l'espressione (15) pud essere riscritta come: 



y-r 



{V 4k*) 



Pirriy-r) 



PPc {<0 x >y\x a >r)*>C PFS X C demod a\--£-;r\p IFT (y-r)e 



(16) 

5 Notiamo che il segnale (16) e significativamente diverso da zero solo nei punti 

aventi y « r , disposti su una retta orizzontale nel dominio (co*; y), cioe parallela alPasse 
delle pulsazioni. Le curve di Fig.9 assumono, dopo questa operazione, 1'aspetto di Fig.10, 
dove sono indicate come Iml, Im2, Im3. In sostanza, Timmagine di ciascun singolo 
riflettore e stata "rimappata", cioe "raddrizzata". 

10 Come si comprende osservaindo la Fig. 10, con Telaborazione sin qui eseguita sul 

segnale generato da un singolo riflettore si e ottenuta, nel dominio delle frequenze, una 
immagine costituita da una linea retta parallela all'asse delle pulsazioni e disposta ad una 
profondita (cioe ad una quota lungo l'asse y) corrispondente alia profondita in cui si trova il 
riflettore. Cioe tutto il contributo generato da un riflettore si trova adesso ad un'unica 

15 profondita, pari a quella dello stesso riflettore. La situazione e equivalente a quella che si 
presentava prima della compressione in profondita. Per ottenere Timmagine dei riflettori e 
a questo punto sufficiente comprimere Timmagine lungo l'asse delle pulsazioni ed eseguire 
Pantitrasformata di Fourier per ricondurre l'immagine nel dominio dello spazio. 

E possibile allora compensare la fase e nello stesso passaggio compensare anche il 

20 fattore introdotto dalla quantita C PFSx C dcmod e dalla antitrasformata di Fourier che 

completera Telaborazione. Grazie alia rimappatura si puo confondere y con r in ogni punto 
deirimmagine e moltiplicare l'espressione (16) per la quantita: 
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co;\ 




4k 2 j 



-1 



(17) 



(18) 



per compensare la fase ed il guadagno, da cui si ottiene il nuovo segnale: 

La fase ed il modulo dipendono linearmente da co x . Questo rende possibile il 
5 calcolo esatto della trasformata inversa di Fourier, da cui si ottiene: 

;pcPc&t>**s}-*ift\-^ (i9) 

dove con airri^l xq , r ) si iridica la trasformata inversa di Fourier della funzione di 
sensibilita laterale a( x; xo , r ). Anch'essa e una funzione a campana, analogamente alia 
funzione di sensibilita laterale, ma molto piu stretta sullo zero. Essa puo essere quindi 

10 considerare una buona approssimazione di una delta di Dirac. 

Dairimmagine descritta dall'espressione (1) o - dall'espressione (2) ad essa 
equivalente, in cui un singolo riflettore dava origine ad una curva di dimensioni 
comparabili con quelle deH'immagine stessa, si e cosi ottenuto una immagine, generata 
dall'espressione (19), che e nulla ovunqiie tranne che in prossimita del punto (x 0 ; r) che 

15 coincide con la posizione del riflettore nello spazio indagato tramite la sonda S. 
L' immagine e stata pertanto correttamente focalizzata. Considerando i tre riflettori che 
inizialmente danno Pimmagine di Fig.8, dopo l'elaborazione di essi si ottiene Timmagine 
schematizzata in Fig.ll. Si noti anche che la fase iniziale <f> Q rimane inalterata attraverso 
tutta l'elaborazione. 

20 L'intera elaborazione e schematizzata e riassunta nel diagramma a blocchi di 

Fig.12. L'immagine ecografica si ottiene eccitando sequenzialmente^uttr^trasduttori che 
formano la cortina di trasduttori della sonda e ricevendo, se(meqfeialmente cbn^ciascun 
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trasduttore, il segnale riflesso. Ricevuti i segnali definenti I'intera immagine ecografica, su 
di essa viene eseguito il procedimento sopra analiticamente descritto, comprendente la 
compressione d'impulso nelle due direzioni eseguita sull'intera immagine. 

L'elaborazione descritta si esegue in questa forma quando il segnale ecografico e 
5 costituito da un segnale di durata relativamente elevata ed a frequenza modulata, ad 
esempio con una modulazione lineare, cioe un cosiddetto segnale chirp. Come risulta dalla 
trattazione sopra riportata, questo consente di ottenere elevate energie del segnale, elevate 
profondita di penetrazione ed una ottima definizione deir immagine. 

Tuttavia, gli ecografi attualmente in commercio utilizzano come segnale di eccitazione 
10 impulsi di breve durata e non modulati in frequenza. In tal caso l'elaborazione risulta molto 
semplificata anche se la qualita deH'immagine che ne deriva e inferiore. Per adattare 
l'algoritmo sin qui descritto a questo tipo di segnali, e sufficiente eliminare i tre passaggi 
della compressione in profondita del segnale (trasformata di Fourier in . profondita, 
compensazione, antitrasformata di Fourier). L'algoritmo assume la forma sinteticamente 
15 rappresentata dal diagramma a blocchi di Fig.13. 

Da un punto di vista matematico, Tespressione del segnale che costituisce il punto 
di partenza in questo caso semplificato non e piu data dall'equazione (2), bensi da 

img(x,t;'x 0 ,r) = a(x;x Q9 r)inv^t-^ 

(20) 

20 Questa espressione si differenzia dall'espressione (2) per l'assenza del termine in cui 
compare la modulazione di frequenza Aft per la diversa espressione della funzione inv(y) r 
che descrive Tinviluppo del segnale ultrasonico trasmesso: in genere questa finestra e 
molto piu stretta del caso precedente ed equivale ad uno o pochi cicli deH'ultrasuono e si 
avvicina essa stessa ad una delta di Dirac. Tutte le restanti funzioni e le costanti coincidono 
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col caso presentato finora. 

Dalla demodulazione coerente e dal passaggio dalla coordinata temporale alia 
coordinata y si ottiene: 

/>/>(x,y;w)-C^ (21) 

Su questo segnale possono essere svolte le stesse considerazioni sul modulo che sono 

state fatte per ilsegnale rappresentato dalFespressione (9), da cui differisce soltanto per la 

funzione p(y) che compare al posto di p JFT (y). 

A partire dall'espressione (21) possono essere ripetuti gli stessi passaggi esposti in 

precedenza ed ottenere alia fine: 

^ / 2& \ ' ■ 

PcP(*> y\ w) = <*ift I - — (* - ^o); r )p(y -r)e- j2kr -e>*» (22) 

Ancora una volta si nota che l'unica differenza rispetto al caso precedente e la 
funzione p(y) al posto di p IFT 

2. Risultati sperimentali: considerazioni generali 

II funzionamento del metodo sopra descritto e stato verificato su acquisizioni reali, 
effettuate sia trasmettendo segnali modulation frequenza, sia con il normale impulso 
trasmesso daU'ecografo. Nel primo caso e stato applicato ralgoritmo completo (Fig. 12), 
nel secondo la versione per segnali non modulati (Fig. 13). In una prima serie di prove. e 
stata ricostruita l'immagine ecografica di un filo metallico teso in acqua, orientato 
ortogonalmente al piano delFimmagine. In una seconda serie di prove e stato utilizzato un 
cosiddetto oggetto campione, che simula riflettori presenti in un reale tessuto da analizzare 
con la tecnica ecografica. Per confronto, le stesse prove sono state eseguite con un ecografo 
funzionante con un sistema di focalizzazione tradizionale. 

Sono stati registrati i segnali provenienti da 250 tracce e 4096 punti per traccia. 
Ogni pixel delle immagini presentate corrisponde a 0,18 mm in senso trasverso e 18,75 \im 
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in profondita; le immagini trattate corrispondono a sezioni larghe 4,5 cm e profonde 7,68 
cm. 

Nelle immagini cui verra fatto riferimento nel seguito la sonda dell'ecografo si 
trova sul lato superiore di ogni immagine ecografica presentata. Come riportato sugli assi, 
5 la coordinata y corrisponde alia profondita, la coordinata x alia direzione trasversa, cioe 
parallela alia lurighezza della cortina dellai sonda. 

Nelle prove eseguite con segnale modulato in frequenza e stato usato un segnale 
ultrasonico di durata 20 >s modulato linearmente in frequenza da 5 MHz a 10 MHz. Per 
generarlo, e stato campionato a 40 MHz e 1 bit un segnale ideale. L'ecografo e stato 

10 riprogrammato in modo da eccitare la sonda con questa sequenza. Grazie alia banda 
limitata dei trasduttori piezoelettrici della sonda, il segnale elettrico con forma d'onda 
quadra diventa un segnale ultrasonico sinusoidale che ben si avvicina a quello ideale. La 
Fig. 14 mostra i segnali utilizzati, per ciascuno dei quali sulla destra della figura e mostrato 
Tingrandimento della porzione evidenziata nel diagramma di sinistra. Piu in particolare in 

15 Fig.l4(a) e mostrato il segnale chirp ideale; in Fig.l4(b) e mostrato il campionamento a 40 
MHz e 1 bit del segnale di Fig. 14(a); in Fig. 14(c) e mostrato il segnale effettivamente 
trasmesso dalla sonda ecografica ed in Fig.l4(d) e rappresentato il segnale ricevuto dalla 
sonda dopo la riflessione su un filo immerso in acqua (l'oggetto utilizzato per la prima serie 
di prove). A causa di riflessioni multiple a diverse profondita, il segnale ricevuto risulta 

20 anche leggermente piu lungo di quello trasmesso. 

II segnale ultrasonico trasmesso mostrato in Fig.l4(c) e stato misurato con un 
idrofono a membrana, tipo 699/1/00002/200 prodotto da Marconi-Caswell, ed un 
oscilloscopio Tektronix TDS 520C. Si puo notare che tale segnale e simile a quello di 
partenza, a parte una modulazione di ampiezza dovuta alia banda non piatta della sonda. II 

25 segnale riflesso e stato ottenuto puntando la sonda verso il filo metallico ed acquisendo il 
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segnale tramite Tecografo. La riflessione degli ultrasuoni ed il cammino del segnale 
elettrico attraverso gli stadi di ricezione aggravant) la modulazione in ampiezza del segnale. 
La modulazione di frequenza rimane invece invariata attraverso tutto il percorso del 
segnale. Come descritto in dettaglio in cio che precede, il primo passo del metodo di 
5 elaborazione demodula coerentemente il segnale utilizzando la frequenza centrale della 
modulazione, nel caso specifico 7,5 MHz. 

Nella serie di prove eseguite con l'impulso normalmente trasmesso dall'ecografo, il 
segnale che eccita la sonda e un singolo ciclo di onda quadra alia frequenza di 6,67 MHz, 
che quindi ha una durata di 150 ns ed e sufficientemente stretto per eccitare la sonda su 

10 tutta la banda disponibile. L'ultrasuono trasmesso segue abbastanza bene l'impulso 
elettrico inviato alia sonda. II segnale ricevuto invece e molto piu lungo del precedente e 
presenta una componente a frequenza piu bassa di quella di eccitazione, circa 4,2 MHz. 
Nella demodulazione coerente della seconda serie di prove e stata usata questa seconda 
frequenza per riportare il segnale in banda base. In Fig.15 sono mostrati i segnali in questo 

15 secondo caso. Piu in particolare, in Fig.l5(a) e rappresentato il segnale elettrico di 
eccitazione della sonda, in Fig.l5(b) il segnale ultrasonico trasmesso, misurato con 
l'idrofono, ed in Fig.l5(c) il segnale ricevuto, acquisito dall'ecografo, costituito dall'eco 
generata dal filo metallico. 

3. Risultati sperimentali: immagine di un filo immerso in acqua 
20 In una prima serie di'prove la sonda ecografica e stata disposta perpendicolarmente 

ad un filo di rame stagnato, con diametro 0,5 mm, teso nell'acqua, in modo che l'ecografia 
ne rilevasse la sezione. II diametro del filo e comparabile con la distanza tra due elementi 
adiacenti della cortina di trasduttori della sonda, che nel caso specificox=e«dLQ,18 mm. Da 



cio si deve prevedere di trovare la sezione del filo contenuta in£un^massim6 Q di<Se tracce 
25 delFecografia, se questa risulta correttamente focalizzata in sensoUiasverso; ^i" 
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Gli ultrasuoni si propagano in acqua alia velocita di circa 1500 m/s, alia frequenza 
di 5 MHz hanno percio una lunghezza d'onda di 0,3 mm, a 10 MHz invece di 0,15 mm. II 
diametro del filo e poco maggiore deila lunghezza d'onda. Questo fa si che soltanto la parte 
superiore del filo rifletta in modo apprezzabile l'onda sonora incidente. Pertanto, il filo, di 
per se gia piccolo, si comporta in modo molto simile ad un riflettore infinitesimo. Ixi 
Fig.l6(a) e mostrato il segnale, registrato dall'ecografo, del filo suddetto. II filo e 
posizionato al centro del pattern e la sonda ecografica e disposta con la propria cortina di 
trasduttori lungo il bordo superiore deirimmagine. Come visibile in Fig. 16(a), il filo genera 
una immagine sfrangiata sia in direzione trasversale, cioe lungo Tasse x, sia nella direzione 
della profondita, cioe in direzione y. In Fig.l6(b) e mostrato un ingrandimento della traccia 
centrale del pattern, passante per il filo. 

La prima operazione deU'elaborazione sul segnale che fornisce 1' immagine di 
Fig. 16(a), la demodulazione coerente, fornisce il risultato di Fig. 17(a), dove e rappresentata 
la parte reale deirimmagine demodulata coerentemente. La Fig.l7(b) riporta 
l'ingrandimento delal traccia centrale del pattern. 

La successiva operazione eseguita suH'immagine demodulata coerentemente e la o 
compressione o focalizzazione in profondita, cioe nella direzione y, operazione il cui 
risultato e visibile in Fig.l8(a) e, in un ingrandimento, in Fig.l8(b), che mostra il modulo 
deirimmagine focalizzata in profondita. Dalle Figg. 18(a) e 18(b) si nota come 
effettivamente la compressione deirimpulso nella direzione y abbia funzionato, generando 
una immagine rappresentata sostanzialmente da un ramo di iperbole con il vertice in 
corrispondenza del punto in cui si trova il filo. La Fig.l9(a) mostra il modulo della traccia 
evidenziata in Figg. 18(a) e 18(b), normalizzato ale la Fig.l9(b) mostra. lo stesso modulo 
normalizzato e presentato in scala logaritmica. Dalla Fig. 19(b) si osserva che la durata del 
segnale compresso in profonditi risulta di circa 15 tracce, considerando 6 dB di dinamica. 
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Eseguita la focalizzazione o compressione in profondita, cioe nella direzione y, 
rimmagine viene sottoposta all'operazione di rimappatura, previa trasformata di Fourier in 
senso trasverso, cioe in direzione x. Le Figg.20(a) e 20(b) mostrano il modulo della 
trasformata di Fourier deU'immagine compressa in profondita, la Fig.20(b) essendo un 
ingrandimento di una porzione della Fig. 20(a). Come visibile in particolare in Fig.20(b), lo 
spettro generato dal filo presenta una marcata curvatura. Rimappando questa immagine 
applicando la trasformazione espressa dalle equazioni (14). Nelle Figg.21(a) e 21(b) (la 
seconda essendo un ingrandimento di una porzione della prima) e riportata 1' immagine 
} ^utteniita/.dispo'ia rimappaturarrora la curvatura dell'irnmagine e siata praticaihence elimmafa 
ed il contributo generato dal filo e riportato ad un'unica profondita, cioe e ricondotto ad 
una immagine sdstanzialmente parallela alia direzione x. L'immagine cosi ottenuta viene 
compressa, cioe focalizzata, in senso trasverso e di essa viene eseguita Pantitrasformata di 
Fourier. II Hsultato di questa ulteriore elaborazione e riprodotta iii.Fig.22(a), che riporta il 
modulo deU'immagine focalizzata. La seziorie del filo viene concentrata in tre tracce, come 
si osserva in particolare nell' ingrandimento riportato nel . riquadro sulla destra 
deU'immagine. La scala dei grigi e lineare. In Fig.22(b) e riportato ancora il modulo 
deU'immagine, con una scala logaritmica dei grigi. Dal nero al bianco sono mostrati i 45 
dB piu significativi. II riquadro a destra e un ingrandimento di una porzione centrales 
deU'immagine, corrispondente alia posizione del filo. Le Figg.23(a) e 23(b) mostrano il 
modulo della traccia orizzdntale, centrata sul filo, deU'immagine focalizzata, normalizzato 
a 1, rispettivamente in scala lineare ed in scala logaritmica. 

Figg.24 e 25 mostrano il segnale reale, acquisito dall'ecografo, e il modulo 
dell'iimmagine focalizzata con una scala di grigi logaritmica, di una immagine analoga 
alia precedente dove, alia sinistra del filo metallico di 0,5 mm di diametro e stato disposto 
un ulteriore filo di nylon con un diametro di 0,1 mm. II filo di nylon, essendo piu piccolo e 
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di materiale diverso, riflette gli ultrasuoni in quantita minore rispetto al filo metallico e 
neirimmagine focalizzata di Fig.25 appare piu piccolo. 

Le Figg.26 a 29(b) mostrano risultati sperimentali ottenuti con lo stesso filo 
utilizzato per le immagini delle figure precedenti, ma con un segnale ecografico 
5 tradizionale, del tipo illustrate in Fig.15. La Fig.26 mostra, analogamente alia Fig. 16, 
Pimmagine o segnale reale acquisito dalPecografo. II filo si trova anche in questo caso al 
centro del pattern. Le Figg.27(a) e 27(b) mostrano, analogamente alle Figg. 23(a) e 23(b) il 
risultato finale dell'elaborazione, cioe il modulo delPimmagine focalizzata, con una scala 
dei grigi lineare e logaritmica rispettivamente. Le Figg.28(a) e 28(b) mostrano, 
10 rispettivamente in scala lineare ed in scala logaritmica, il modulo della traccia orizzontale, 
centrata sul filo, delPimmagine focalizzata, normalizzato a 1. Tali figure equivalgono alle 
Figg.23(a) e 23(b) ottenute con il segnale ecografico modulato in frequenza. Infine, le 
Figg.29(a) e 29(b) mostrano, rispettivamente in scala lineare ed in scala logaritmica, il 
modulo normalizzato a 1 della traccia verticale centrata sul filo delPimmagine focalizzata. 
15 Tali figure equivalgono alle Figg. 19(a) e 19(b) precedenti. 

Si osserva che il metodo di elaborazione della presente invenzione consente di 
ottenere la desiderata focalizzazione trasversa delPimmagine, anche in assenza di 
modulazione di frequenza nel segnale ecografico di partenza. I risultati sono, infatti, simili 
a quelli conseguiti con un segnale modulato in frequenza (segnale chirp), bench e si noti un 
20 peggioramento della risoluzione in profondita. 

Le Figg. 30(a) e 30(b) mostrano, rispettivamente in scala lineare di grigi ed in scala 
logaritmica di grigi, il risultato di una focalizzazione tradizionale delPimmagine dello 
stesso oggetto. Rispetto ai risultati ottenuti con il metodo della presente invenzione ed 
/ v " • esposti rielle figure precedenti, il filo emerge maggiormente dal fondo. Si ottiene, cioe, 
. 25 un^immagine a maggiore contrasto. Tuttavia, la risoluzione ottenuta e peggiore ed inoltre 
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risulta non uniforme. Cio si puo osservare con particolare evidenza nelle Figg.31(a), 31(b) 
e 31(c) e nelle Figg.32(a), 32(b) e 32(c). 

La Fig.31(a) mostra l'andamento del modulo della traccia orizzontale centrata sul 
filo, normalizzato a 1 ed in scala logaritmica ottenuta dalFecografo tradizionale. La 
5 Fig.31(b) mostra un ingrandimento (curva a) della porzione nel riquadro di Fig.31(a) 
sovrapposta all'andamento corrispondente ottenuto con il metodo dell' invenzione 
utilizzando un segnale non modulato in frequenza (curva b). La Fig.31(c) mostra la stessa 
curva a di Fig.31(b) confrontata con l'andamento corrispondente (curva b) ottenuto con il 
metodo della presente invenzione e con segnale ecografico modulato in frequenza (segnale 
10 chirp). 

Le Figg.32(a), 32(b), 32(c) mostrano un analogo confronto dell'andamento dei 
moduli delle tracce verticali centrate sul filo. In Fig.32(a) e mostrato l'andamento nel caso 
di ecografo tradizionale; in Fig. 32(b) tale andamento e confrontato con quello ottenuto con 
il metodo di focalizzazione della presente invenzione utilizzando un segnale ecografico 
15 tradizionale, cioe non modulato in frequenza. La Fig.32(c), infine, mostra il confronto con 
il modulo della corrispondente traccia verticale ottenuta con un segnale chirp, cioe 
modulato in frequenza. 
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Rivendicazioni 

1. Un metodo per la focalizzazione di immagini ecografiche, comprendente le 

fasi di: 

• inviare in un volume sotto indagine una serie di segnali ultrasonici di eccitazione, 
tramite una cortina di trasduttori allineati lungo una direzione trasversale (x), detti 
segnali ultrasonici propagandosi in profondita in detto volume secondo una 
direzione di propagazione (y); 

• raccogliere, tramite detti trasduttori, segnali riflessi da riflettori disposti nel volume 
sotto indagine; 

• eseguire su detti segnali riflessi una trasformata in direzione trasversale da un 
dominio spaziale (x,y)> definito da detta direzione trasversale (x) e da detta 
direzione di propagazione (y), ad un primo dominio trasformato; 

• applicare, nel dominio trasformato, una trasformazione bidimensionale, per 
raddrizzare ogni immagine curva (Ipl, Ip2, Ip3) di un riflettore in detto volume 
sotto indagine e renderla sostanzialmente ortogonale alia direzione di propagazione 

(y); • 

• comprimere, nella direzione trasversale (x) ciasciina curva raddrizzata (Iml, Im2, 
Im3) per concentrare detta immagine raddrizzata in. una zona centrata in 
corrispondenza della posizione di detto riflettore lungo detta direzione trasversale 

(*)• 

2. Metodo come da rivendicazione 1, in cui sui segnali raccolti viene eseguita 
una demodulazione coerente del segnale per riportarlo in banda base. ; 

3. Metodo come da rivendicazione 1 o 2, in cui detta trasformata e un 
trasformata di Fourier in direzione trasversale (x). 

><4. Metodo come da rivendicazione 1 o 2 o 3, in cui la compressione trasversale 
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di dette curve raddrizzate viene completata con una compensazione del guadagno e di fase 
e con una antitrasformata, che riporta l'immagine dal dominio della frequenza al dominio 
dello spazio. 

5. Metodo come da rivendicazione 1, 2. 3 o 4 in cui detta trasformazione 
bidimensionale e espressa da: 

■ (14) 
dove co x e la coordinata nella direzione delle frequenze 

k e la costante di propagazipne pari a ^> 

' • A, , 

e X e la lunghezza d'onda del segnale ultrasonico trasmesso. 

6. Metodo come da una o piii delle rivendicazioni precedent^ Jn cui detti 
segnali ultrasonici di eccitazione sono modulati in frequenza. 

7. Metodo come da rivendicazione 6, in cui detta modulazione in frequenza e 
t una modulazione lineare. 

8. . Metodo come da rivendicazione 6 o 7, in cui detti segnali ultrasonici di 
eccitazione presentano un inviluppo rettangolare. 

9. Metodo come da rivendicazione 6 o 7 od 8, in cui i segnali riflessi da detti 
riflettori vengono riportati in banda base tramite una demodulazione coerente. 

10. Metodo come da rivendicazione 6 o 7 od 8 o 9, in cui i segnali riflessi da 
detti riflettori nel volume sotto indagine vengono compressi nella direzione di 
propagazione (y). 

11. Metodo come da rivendicazione 10, in cui detti segnali riflessi vengono 
compressi nella direzione di propagazione tramite una trasformata nella direzione di 

s 

- v » 

Pag. 35 di 36 

■C:\documenti_michele\descrizioni\actl91102.doc 



Ufficio Tecnic Ing. A. Mannucci srl 

propagazione (y), dal dominio spaziale ad un dominio delle frequenze (r, co y ) y una 
compensazione della fase e del guadagno, ed una successiva trasformata inversa. 

12. Metodo come da rivendicazione 11, in cui delta trasformata nella direzione 
di propagazione e una trasformata di Fourier. 
5 13. Metodo come da una o piu delle rivendicazioni precedent^ in cui detti 

segnali ultrasonici di eccitazione vengono inviati in detto volume sotto indagine 
sequenzialmente da singoli trasduttori o da gruppi comprendenti un numero limitato di 
trasduttori, ed in cui Pimmagine del volume sotto indagine viene ottenuta raccogliendo 
sequenzialmente, per ciascun segnale inviato da ciascun singolp trasduttore o da ciascun 

10 gruppo di trasduttori, i segnali riflessi da detti riflettori tramite e raccolti dal trasduttore o 
dal gruppo di trasduttori che ha emesso il relativo segnale inviato. 

14. Un sistema ecografico comprendente almerio una sonda con una cortina di 
trasduttori allineati secondo una direzione trasversale di allineamento, mezzi per eccitare 
sequenzialmente detti trasduttori, mezzi di elaborazione per ricevere ed elaborare segnali 

15 riflessi da riflettori contenuti in un volume sotto indagine in cui vengono propagati segnali 
ultrasonici generati da detti trasduttori, detti -mezzi di elaborazione eseguendo la 
focalizzazione in direzione trasversale del segnale ecografico ricevuto da detti trasduttori 
con un metodo secondo una o piu delle rivendicazioni precedents 

FIRENZE 2 1 MOV, 2002 
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